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We  present  a  molecular-level  simulation  study  of  the  effects  of  confinement  on  chemical  reaction 
equilibrium  in  nanoporous  materials.  We  use  the  reaction  ensemble  Monte  Carlo  (RxMC)  method  to 
investigate  the  effects  of  temperature,  nanopore  size,  bulk  pressure,  and  capillary  condensation  on 
the  nitric  oxide  dimerization  reaction  in  a  model  carbon  slit  nanopore  in  equilibrium  with  a  bulk 
reservoir.  In  addition  to  the  RxMC  simulations,  we  also  utilize  the  molecular-dynamics  method  to 
determine  self-diffusion  coefficients  for  confined  nonreactive  mixtures  of  nitric  oxide  monomers 
and  dimers  at  compositions  obtained  from  the  RxMC  simulations.  We  analyze  the  effects  of  the 
temperature,  nanopore  width,  bulk  pressure,  and  capillary  condensation  on  the  reaction  equilibrium 
with  respect  to  the  reaction  conversion,  fluid  structure,  and  self-diffusion  coefficients.  We  show  that 
the  influence  of  the  temperature,  nanopore  size,  and  capillary  condensation  on  the  confined  reaction 
equilibrium  is  quite  dramatic  while  the  effect  of  the  bulk  pressure  on  the  reaction  equilibrium  in  the 
carbon  slit  nanopore  is  only  moderate.  This  work  is  an  extension  of  previous  work  by  Turner  et  al. 
[J.  Chem.  Phys.  114,  1851  (2001)]  on  the  confined  reactive  nitric  oxide  system. 
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I.  INTRODUCTION 

With  the  rapid  growth  of  nanotechnology  and  the  inven¬ 
tion  of  various  nanoporous  materials,  some  of  these  materials 
have  been  proposed  as  vehicles  for  nanochemical  devices 
such  as  nanoscale  reactors  and  nanoscale  membrane 
reactors.1  However,  further  development  of  these  applica¬ 
tions  is  impossible  without  a  fundamental  understanding  of 
the  thermodynamic,  reaction,  adsorption,  and  transport 
mechanisms  of  fluids  confined  in  nanopores. 

We  now  have  a  fairly  good  understanding  of  the  influ¬ 
ence  of  confinement  on  the  physical  properties  of 
nanophases  for  simple  and  moderately  complex  nonreactive 
fluids."'  It  is  well  established,  for  example,  that  confinement 
brings  about  drastic  changes  in  the  thermodynamic  proper¬ 
ties  of  nonreactive  fluids  such  as  narrowing  the  vapor-liquid 
coexistence  curve,  lowering  the  pore  critical  temperature,  in¬ 
creasing  the  average  fluid  densities  in  pores,  and  the  appear¬ 
ance  of  new  types  of  phase  transitions  not  found  in  the  bulk 
phase. 
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By  contrast,  significantly  less  is  known  about  the  effects 
of  confinement  on  chemical  properties,  particularly  on 
chemical  reaction  equilibria.  A  chemical  reaction,  confined  to 
a  nanosized  environment,  can  have  a  different  outcome  as 
compared  to  the  same  reaction  in  the  bulk  phase.  For  ex¬ 
ample,  generally  the  nanopore  phase  has  a  higher  density 
than  the  corresponding  bulk  phase;  based  upon  Le  Chat- 
elier’s  principle  this  results  in  an  increase  in  yield  for  reac¬ 
tions  in  which  there  is  a  decrease  in  the  total  number  of 
moles.  Conversely,  a  drop  in  yield  occurs  in  reactions  where 
the  number  of  moles  increases.  Further,  some  components  of 
the  reactive  mixture  are  selectively  adsorbed  on  the  solid 
surfaces,  subsequently  affecting  the  reaction  equilibrium. 
Still  further,  molecular  orientations  can  be  strongly  influ¬ 
enced  by  proximity  to  a  solid  surface  which  can  also  shift  the 
reaction  equilibrium  relative  to  the  bulk  phase  equilibrium. 
Finally,  phase  transitions  such  as  capillary  condensation  are 
expected  to  have  a  strong  influence  on  the  reaction  conver¬ 
sion  in  the  nanopores. 

Molecular-level  simulation  studies  of  the  effects  of  con¬ 
finement  on  reaction  equilibria  were  pioneered  by  Borowko 
et  al.4  and  Borowko  and  Zagorski5  for  model,  reversible  re¬ 
actions  in  slit  pores  and  by  Turner  et  a/.6-8  for  realistic,  re¬ 
versible  reactions  (nitric  oxide  dimerization  and  ammonia 
synthesis)  in  carbon  micropores  and  carbon  nanotubes.  The 
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study  by  Turner  et  al.6  was  inspired  by  the  experimental 
work  of  Kaneko  et  al.9  and  Nishi  et  al.10  who  attempted  to 
experimentally  measure  reaction  equilibria  in  carbon  nanop¬ 
ores.  More  recently,  Peng  et  al. 11  simulated  chemical  reac¬ 
tion  equilibria  of  ammonia  synthesis  in  two  porous  materials 
(MCM-41  and  pillared  clays)  and  Hansen  et  al.x~  simulated 
the  reaction  equilibria  of  a  metathesis  reaction  system  in 
silicalite-1  pores. 

All  the  above-mentioned  simulations  mainly  focused  on 
the  influence  of  the  pore  size  on  the  reaction  conversion.  In 
addition,  they  also  considered  the  effects  of  chemical  and 
physical  surface  heterogeneities  on  the  reaction  conversion, 
as  well  as  the  influence  of  the  variation  of  the  equilibrium 
constant  across  the  pore  on  the  reaction  equilibrium.  These 
authors  employed  the  reaction  ensemble  Monte  Carlo 
(RxMC)  simulation  technique13-15  that  enables  one  to  di¬ 
rectly  simulate  the  equilibrium  properties  of  chemically  re¬ 
acting  systems.  The  method  requires  only  a  knowledge  of  the 
reaction  species  intermolecular  potentials  and  their  ideal-gas 
properties,  in  addition  to  specifying  the  system  stoichiometry 
and  the  thermodynamic  constraints. 

While  dramatically  influencing  the  thermodynamic  prop¬ 
erties,  confinement  also  affects  the  transport  properties  of 
fluids  inside  nanopores.  Moreover,  confinement  contributes 
significantly  to  both  the  thermodynamic  and  transport  prop¬ 
erties  if  chemical  reactions  and  mixture  separation  are  occur¬ 
ring  simultaneously  in  the  nanoporous  materials.  Physical 
space  restrictions  based  on  the  fluid  particle  size  or  structure 
may  limit  the  flow  of  particles  through  particular  nanopores 
in  the  material.  Furthermore,  attraction  of  the  nanopore  sur¬ 
face,  i.e.,  physisorption,  may  play  a  critical  role.  Even  further 
complicating  matters,  fluid  particles  may  chemisorb.  In  ad¬ 
dition  to  phenomena  occurring  between  fluid  particles  and 
the  nanopore  surface,  the  transport  behavior  in  the  nano¬ 
porous  material  becomes  increasingly  more  complex  if 
chemical  reactions  occur  between  fluid  particles. 

Recently,  we  proposed  a  novel  simulation  tool  for  such 
scenarios,  termed  the  dual  control  cell  reaction  ensemble 
molecular-dynamics  (DCC-RxMD)  method.16'17  The  DCC- 
RxMD  method  simulates  fluid  mixtures  that  are  simulta¬ 
neously  chemically  reacting  and  adsorbing  in  a  nanoporous 
material  by  coupling  a  nonequilibrium  molecular-dynamics 
method  with  the  RxMC  and  grand  canonical  ensemble 
Monte  Carlo  methods.  The  simulation  setup  mimics  an  ex¬ 
perimental  system  where  the  control  cells  that  maintain  the 
desired  reaction  and  flow  conditions  are  in  direct  physical 
contact  with  the  adsorbing  media.  The  DCC-RxMD  method 
was  applied  to  the  dry  reforming  of  methane  reaction  in 
model  nanoscale  membrane  reactors.16  These  DCC-RxMD 
simulation  studies  provided  unique  insight  into  the  effects  of 
the  membrane  structure  and  porosity  and  imposed  pressure 
gradients  on  the  reaction  conversion. 

In  addition  to  these  simulation  and  experimental  studies, 
density-functional  theory  (DFT)  studies  on  model  dimeriza¬ 
tion  reactions  in  slit  pores  have  also  been  performed.13  Trip- 
athi  and  Chapman  extended  the  DFT  method  developed  by 
Chapman  and  co-workers  for  confined  associating 
systems19'20  to  the  study  of  confined  dimerization  reactions. 
They  studied  the  effects  of  the  pore  size  and  capillary  con¬ 


densation  on  the  reaction  equilibrium.  DFT  predicts  an  in¬ 
crease  in  the  reaction  conversion  due  to  capillary  condensa¬ 
tion,  with  the  impact  being  most  significant  in  the  smaller 
pores. 

In  this  paper,  we  continue  the  work  of  Turner  et  al.6  on 
the  nitric  oxide  dimerization  reaction,  2NO=^(NO)2,  in 
model  carbon  slit  nanopores.  However,  motivated  by  the 
findings  of  the  DFT  studies  by  Tripathi  and  Chapman,18  we 
shift  our  focus  to  the  effects  of  bulk  pressure  and  capillary 
condensation  on  the  confined  reaction  equilibrium  and  fluid 
structure.  To  provide  further  insight,  we  report  self-diffusion 
coefficients  obtained  from  MD  simulations  of  confined  non¬ 
reactive  mixtures  of  nitric  oxide  monomers  and  dimers  at 
compositions  determined  from  RxMC  simulations.  More  in- 
depth  studies  of  the  transport  properties  of  this  system  will 
be  considered  in  subsequent  work.  For  these  future  studies 
the  implementation  of  the  aforementioned  DCC-RxMD 
method  will  be  necessary. 

Besides  being  studied  previously,  the  nitric  oxide  dimer¬ 
ization  reaction  is  interesting  for  a  number  of  reasons.  The 
reaction  is  an  exothermic,  thermodynamically  driven  reac¬ 
tion  in  which  there  is  a  decrease  in  the  total  number  of 
moles.  The  reaction  is  important  in  atmospheric  chemistry  as 
well  as  in  the  human  body  where  it  regulates  blood  pressure. 
Moreover,  predicting  the  effects  of  confinement  on  NO 
dimerization  is  critical  to  pollution  abatement  since  activated 
carbons  are  commonly  used  for  the  removal  of  nitrogen  ox¬ 
ides  from  autoexhaust  and  industrial  effluent  gas  streams. 
From  a  modeling  standpoint,  the  NO  dimerization  reaction  in 
carbon  slit  nanopores  is  an  ideal  reaction  to  simulate.  Carbon 
nanopores  exhibit  quite  strong  solid-fluid  interactions  ’  and 
the  molecules  involved  in  the  reaction,  NO  and  (NO)2,  are 
simple  and  easy  to  model.  Due  to  the  NO  paramagnetism  and 
(NO)2  diamagnetism,  the  composition  of  the  reactive  mix¬ 
ture  in  activated  carbons  can  be  measured  by  magnetic 
susceptibility,9  providing  experimental  data  for  comparison. 

The  paper  is  organized  as  follows.  The  simulation  meth¬ 
odology  is  described  in  Sec.  II  along  with  the  simulation 
details.  Molecular  models  for  fluid-fluid  and  solid-fluid  inter¬ 
actions  together  with  the  evaluation  of  the  equilibrium  con¬ 
stant  are  presented  in  Sec.  III.  Results  are  discussed  in  Sec. 
IV.  Finally,  Sec.  V  gives  our  conclusions. 


II.  SIMULATION  METHODOLOGY 

In  this  work,  we  used  the  RxMC  method  to  simulate  the 
reaction  equilibrium  of  the  2NO=^(NO)2  system  in  a  model 
carbon  slit  nanopore  in  equilibrium  with  a  bulk  reservoir. 
Constant-temperature-constant-volume  ( NVT )  MD  simula¬ 
tions  were  carried  out  to  determine  the  self-diffusion  coeffi¬ 
cients  of  the  corresponding  confined  NO/(NO)2  equilibrium 
mixture.  We  assessed  the  dependency  of  the  reaction  equilib¬ 
ria  and  self-diffusion  coefficients  on  various  factors,  includ¬ 
ing  the  nanopore  width,  bulk  reservoir  pressure,  temperature, 
and  capillary  condensation.  For  each  case,  the  fluid  structure 
was  analyzed  by  constructing  density  profiles  of  each  species 
in  the  nanopore. 
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A.  Reaction  ensemble  Monte  Carlo 

13—15  r 

The  RxMC  method  is  a  powerful  simulation  tool  for 
studying  chemically  reacting  mixtures.  The  method  only  re¬ 
quires  inputting  the  intermolecular  potentials  and  the  ideal- 
gas  properties  of  the  reaction  species  that  are  present.  Most 
notably,  the  method  does  not  require  a  reactive-type  potential 
that  mimics  bond  breakage  and  formation,  see  e.g.,  Refs.  21 
and  22.  Reactions  are  simulated  by  performing  forward  and 
reverse  reaction  steps  according  to  the  RxMC  algorithm 
which  guarantees  that  the  reaction  equilibrium  criteria  for  a 
set  of  R  linearly,  independent  chemical  reactions 

C 

2vb  =  °,  7  =  1,2,  ...,R  (1) 

i—  1 

are  established."3  In  Eq.  (1),  Vjt  is  the  stoichiometric  coeffi¬ 
cient  of  species  i  in  chemical  reaction  j  and  /x,  is  its  chemical 
potential.  The  reaction  equilibrium  condition  for  our  NO 
dimerization  reaction  2 NO  ^  (NO)2  is  then 

AhNoq  -  2/xn0  =  0 .  (2) 


Forward  and  reverse  reaction  steps  are  accepted  with  prob¬ 
abilities 


min 


,  T  ANO(/VNO  -  1) 

1 , - expl 

y  ^(NO)2+l 


A  U\ 

kBT) 


(3) 


and 


spherical  cutoff  radius  equal  to  half  the  box  length  was  used 
and  the  long-range  correction  for  the  configurational  energy 
was  included,  assuming  that  the  radial  distribution  function 
is  unity  beyond  the  cut-off  radius. 

The  simulations  were  initiated  by  randomly  placing  NO 
molecules  into  the  nanopore  and  bulk  simulation  boxes.  The 
initial  number  of  NO  molecules  was  chosen  such  that  a  sta¬ 
tistically  reasonable  number  of  molecules  (—500)  were 
present  in  the  bulk  phase  once  the  system  had  equilibrated. 
The  RxMC  simulations  were  organized  in  cycles  as  follows. 
Each  cycle  consisted  of  four  steps:  nD  displacement  moves 
(translation  or  rotation),  nv  volume  moves  for  the  bulk  phase, 
rij:  reaction  moves,  and  nT  NO  molecule  transfers.  The  four 
types  of  moves  were  selected  at  random  with  fixed  probabili¬ 
ties,  chosen  such  that  the  ratio  nD:nv:n^:nT  in  each  cycle 
was  TV :  1 :7V:  TV,  where  N  was  about  10%-20%  greater  than 
the  equilibrium  number  of  molecules  during  a  simulation 
run.  The  acceptance  ratios  for  translational  and  rotational 
moves,  and  for  volume  changes,  were  adjusted  to  be  approxi¬ 
mately  30%.  After  an  equilibration  period,  we  typically  gen¬ 
erated  6X  104  cycles  to  accumulate  averages  of  the  desired 
quantities.  The  precision  of  the  simulated  data  was  obtained 
using  block  averages,  with  5000  cycles/block.  In  addition  to 
ensemble  averages  of  the  quantities  of  direct  interest,  we  also 
carefully  monitored  the  convergence  profiles  of  the  thermo¬ 
dynamic  quantities  as  the  system  traveled  through  phase 
27 

space. 


min 


V  N(NQ)2  j 

r  (Nn  o  +  i)(a?No  +  2) 


A U\ 

wJ’ 


(4) 


respectively.  In  Eqs.  (3)  and  (4),  V  is  the  system  volume,  T  is 
the  temperature,  A,-  is  the  number  of  particles  of  species  i,  kB 
is  Boltzmann’s  constant,  1  is  the  ideal-gas  quantity  defined 
as 


r  = 


kBT 

-\k, 

po 


(5) 


P°  is  the  standard-state  pressure  (taken  to  be  1  bar),  AU  is 
the  change  in  the  configurational  energy  U  due  to  forward 
and  reverse  reaction  attempts,  and  K  is  the  equilibrium  con¬ 
stant. 

Reaction  steps  and  standard  Monte  Carlo  displacement 
moves"  '  were  carried  out  in  both  the  bulk  and  nanopore 
phases.  The  simulation  of  a  nanopore  phase  in  equilibrium 
with  a  bulk  reservoir  phase  requires  volume  changes  for  the 
bulk  phase  along  with  particle  interchanges  between  the  na¬ 
nopore  and  bulk  phases.26  Particle  exchanges  between  the 
phases  were  only  performed  for  the  NO  monomer,  which  is 
the  smaller  of  the  molecules  in  the  dimerization  reaction.6 

The  nanopore  consisted  of  two  structureless  confining 
walls  separated  by  a  distance  H  in  the  z,  direction  with  peri¬ 
odic  boundary  conditions  applied  in  the  x  and  y  directions 
only.  The  bulk  phase  was  represented  by  a  cubic  simulation 
box  where  the  minimum  image  convention  and  periodic 
boundary  conditions  were  applied.  In  the  case  of  the  nanop¬ 
ore  phase,  a  cutoff  equal  to  half  the  maximum  box  size  was 
utilized  where  the  long-range  correction  for  the  configura¬ 
tional  energy  was  ignored.  In  the  case  of  the  bulk  phase,  a 


B.  Molecular  dynamics 

After  determining  the  reaction  equilibrium  concentra¬ 
tions  in  the  nanopore  phase  from  RxMC  simulations,  the 
transport  properties  of  these  fluid  mixtures  were  determined 
using  the  following  approach.  A  fluid  mixture  matching  the 
final  average  concentrations  from  the  RxMC  simulation  was 
created  in  the  nanopore.  We  then  employed  a  standard  NVT 
MD  method"4  "3  to  simulate  particle  movements  trough  time 
space.  Since  a  smooth  potential  was  used  to  represent  the 
graphite  planes  (as  opposed  to  an  atomistically  detailed  rep¬ 
resentation),  we  needed  to  account  for  the  exchange  of  mo¬ 
mentum  which  would  take  place  between  fluid  molecules 
and  the  carbon  atoms.  Therefore,  within  the  MD  method 
framework  we  implemented  the  so-called  diffuse  boundary 
condition  based  on  the  work  of  Travis  and  Gubbins  as  fol¬ 
lows.  After  each  MD  time  step  we  monitored  whether  the 
following  two  conditions  were  satisfied:  (i)  the  center-of- 
mass  momentum  of  a  given  molecule  in  the  direction  perpen¬ 
dicular  to  the  wall  has  reversed  in  sign  and  (ii)  the  center-of- 
mass  momentum  of  that  same  molecule  was  within  the 
repulsive  region  of  the  solid-fluid  potential.  If,  and  only  if, 
both  of  these  conditions  were  satisfied,  then  we  reselect  the 
center-of-mass  momentum  of  that  molecule  in  the  directions 
parallel  to  the  confining  walls  from  a  Maxwell-Boltzmann 
distribution  corresponding  to  the  prescribed  temperature."4 
Such  an  approach  does  not  alter  the  angular  velocity  of  mol¬ 
ecules  that  are  approximately  spherical-like  and  diffusing  at 
moderate  to  high  temperatures  (both  approximations  are  sat¬ 
isfied  in  this  study). 

The  equations  of  motion  in  the  NVT  MD  method  were 
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solved  using  the  leap-frog  algorithm'  and  by  implementing 
the  damped  force  method  of  Brown  and  Clark  and 
Fincham30  along  with  the  Berendsen  thermostat31  to  maintain 
constant  temperature.  A  time  step  of  At  =  2.5  fs  and  a  cou¬ 
pling  time  constant  77-=  0.01  ps  in  the  Berendsen  thermostat 
were  used  for  a  total  simulation  time  of  1  ns  following  an 
equilibration  period. 

During  the  course  of  the  NVT  MD  simulation,  we  evalu¬ 
ated  mean-square  displacements  in  directions  parallel  to  the 
walls,  denoted  as  MSD^,  where  a=  (x,y)  and  i 
=  [NO,(NO)2].  MSD'a  values  were  used  to  calculate  self¬ 
diffusion  coefficients  in  one  dimension,  D'a,  by  means  of  the 
Einstein  equation,24  where 


D'  =  lim 


MSD'a 
2 t 


and 


MSD^=(Aa'(f)> 


1  N‘ 


0=1 


a){t) 


■  «}(0)|2 


(6) 


(7) 


In  Eqs.  (6)  and  (7),  t  is  the  time,  cr'  is  a  component  of  the 
unfolded  position  of  a  molecule  j  (i.e.,  not  subject  to  periodic 
boundary  conditions),  and  (■)  is  the  appropriate  ensemble 
average.  In  a  slit  nanopore,  D'x  and  D‘  should  be  the  same 
within  simulation  uncertainties  and  we  thus  report  one¬ 
dimensional  self-diffusion  coefficients,  D\  as  average  values 
of  the  corresponding  D\  and  D\  quantities. 


III.  MOLECULAR  MODEL 
A.  Fluid-fluid  interaction 

As  in  the  previous  simulation  studies  of  the  NO  dimer¬ 
ization  reaction,6'14  we  used  the  model  proposed  by  Kohler  et 
al.  to  describe  intermolecular  interactions  in  the  mixture  of 
NO  and  (NO)2.  The  model  treats  NO  as  a  single  Lennard- 
Jones  (LJ)  sphere  and  (NO)2  as  a  two-site  LJ  molecule  with 
the  bond  length  /  equal  to  the  experimental  value  of 
0.2237  nm.  For  the  monomer,  the  model  uses  the  LJ  energy 
parameter  s/kB=  125.0  K  and  the  LJ  size  parameter  a 
=  0.317  15  nm.  The  individual  LJ  parameters  for  each  site  in 
the  dimer  are  the  same  as  those  for  the  monomer.  Due  to  the 
weak  dipole  and  quadrupole  moments,  the  model  neglects 
electrostatic  forces. 


B.  Solid-fluid  interaction 

When  simulating  equilibrium  fluid  properties  in  carbon 
slit  nanopores,  it  is  possible  to  treat  the  walls  as 
structureless.  In  this  case,  a  solid-fluid  intermolecular  poten¬ 
tial  is  obtained  by  replacing  the  sum  over  solid-fluid  particle 
interactions  by  a  sum  of  integrals  over  wall  atoms  in  a  given 
plane.  This  is  a  reasonable  approximation  when  the  fluid 
molecule  is  large  compared  with  the  spacing  between  wall 
atoms.  In  graphitic  carbons,  the  C-C  spacing  between  sur¬ 
face  carbon  atoms  is  about  0.142  nm  and  molecules  such  as 
NO  (cr=0.317  15  nm)  or  (NO)2  (cr=0.317  15  nm,/ 

=  0.2237  nm)  feel  only  a  rather  small  corrugation  in  the 
solid-fluid  interaction  as  it  moves  parallel  to  the  surface.  As¬ 
suming  a  LJ  potential  for  the  wall  atom-fluid  interaction  and 


-  H  =  0.8 

nm 

- H  =  1.0 

nm 

- H  =  1.4 

nm 

.  H  =  1.7 

nm 

- -  H  =  2.0 

nm 

FIG.  1.  Solid-fluid  interaction  potential,  us{,  for  the  NO  monomer  in  the 
model  carbon  slit  nanopore  of  various  widths  H;  z  is  the  distance  of  the  NO 
molecule  from  the  center  of  the  nanopore  which  is  located  at  z=0  and  kB  is 
Boltzmann’s  constant. 


integrating  over  the  interactions  with  individual  carbon  at- 
oms  in  each  graphite  plane,  Steele  obtained  the  10-4-3  po¬ 
tential  for  the  interaction  of  a  fluid  atom  with  the  graphite 
wall, 

«sffe)  =  27resf°2S!P.A 

(8) 


2/gsfV0  /  4 

5  \z)  \z)  3A(z  +  0.61A)3 


In  Eq.  (8),  7  is  the  distance  of  the  fluid  atom  from  the  graph¬ 
ite  surface,  esf  and  crsf  are  the  LJ  potential  parameters  for  the 
solid-fluid  interaction,  A  =  0.335  nm  is  the  interplanar  spac¬ 
ing  in  graphite,  and  ps=  114  nm”3  is  the  number  density  of 
carbon  atoms.  The  ssf  and  <xsf  parameters  are  obtained  from 
the  Lorentz-Berthelot  mixing  rules24  with  the  carbon  LJ  po¬ 
tential  parameters  ess/kB  =  28  K  and  crss=0.340  nm. 

For  a  slit  nanopore  of  width  H,  the  fluid  molecule  inter¬ 
acts  with  both  graphite  walls  so  that  the  total  solid-fluid  in¬ 
teraction  is  given  as  the  sum  of  us[(z)  and  usf(H-z).  As  an 
example,  the  solid-fluid  interaction  potential  for  the  NO 
monomer  in  the  model  carbon  slit  nanopore  of  various 
widths  is  shown  in  Fig.  1. 


C.  Equilibrium  constant 


The  RxMC  method  requires  knowledge  of  the  equilib¬ 
rium  constant  K.  For  the  NO  dimerization  reaction 
2NO; 


j23 


(NO)2,  K  is  defined” 

-  2/4 Iio  \ 
RT  j" 


K  =  expl  - 


0 

/Ano), 


(9) 


where  /j^(T,P°)  is  the  molar  standard  chemical  potential  of 
species  i  at  T  and  P°,  and  R  is  the  universal  gas  constant.  The 
dependence  of  /j^(T,P°)  on  P°  is  usually  suppressed,  and 
p%T)  may  be  expressed  as 
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Jll%T)=Ii6(T)-Ts%T), 

(10) 

where  the  molar  enthalpy,  h®(T),  and 
s6(T),  may  be  expressed  as 

the  molar  entropy, 

h?(T)=AHfi(Tr)  +  f  c°pi(T)dT, 

Tr 

(ID 

°lT\  °(T  1  j-  f  Cpi ^  IT 

s,  (T)  -  s,  (Tr)  +  1  dT. 

Tr 

(12) 

In  Eqs.  (11)  and  (12),  A Hfi(Tr)  is  the  enthalpy  of  formation 
of  species  i  at  the  reference  temperature  Tr  (taken  to  be 
298.15  K)  and  cpi  is  its  ideal-gas  heat  capacity. 

In  general,  K  refers  to  the  isolated-molecule  (ideal-gas) 
properties.  The  most  accurate  values  for  these  properties  are 
available  in  thermochemical  tables  such  as  the  TRC  or 
JANAF  tables,34’35  which  are  constructed  using  the  most  ac¬ 
curate  available  partition  function  (PF)  data.  Hence  in  this 
work,  values  of  K  at  different  7”s  were  evaluated  by  means 
of  the  TRC  thermochemical  table34  and  Eqs.  (9)— ( 12).  In 
contrast.  Turner  et  al.6  and  Johnson  et  a/.14  in  their  RxMC 
study  of  the  NO  dimerization  system  calculated  K  using  the 
internal  degrees  of  freedom  of  the  NO  and  (NO)2  molecules 
(translational,  rotational,  vibrational,  and  electronic)  based 
on  classical  statistical  mechanics.36  In  principal,  both  ap¬ 
proaches  should  yield  the  same  K  values  since  fi-  is  related 
to  the  internal  PF,  (Ref.  37)  via 


fj°  =  -RT  In' 


r  v  e  |  rrX 

kBT\ 
A3  P°j 


(13) 


In  Eq.  (13),  A,-  is  the  de  Broglie  thermal  wavelength,  and  qr{, 
q ",  and  q'\  are,  respectively,  rotational,  vibrational,  and  elec¬ 
tronic  contributions  to  the  internal  PF,.  However,  due  to  a 
large  variation  in  the  literature  values  of  the  (NO)2  dissocia¬ 
tion  energy  D0,  which  is  needed  for  the  evaluation  of  ^(N0>  , 
Turner  et  al.6  and  Johnson  et  al.U  adjusted  the  value  of  D0 
from  a  series  of  short  simulation  runs  performed  at  110  and 
116  K,  respectively.  D0  was  adjusted  until  a  reasonable 
agreement  with  the  experimental  data  of  dimer  concentra¬ 
tions  in  the  bulk  saturated  liquid  was  found.  Turner  et  al. 
obtained  Z)0=  10.9  kJ/mol  at  T=  I  I  O  K  and  Johnson  et  al. 
obtained  D0=  13.6  kJ/mol  at  T=116K;  both  values  are 
within  the  wide  range  of  D0  values,  7.64-13.8  kJ/mol,  re¬ 
ported  in  the  literature.39  42 


IV.  RESULTS  AND  DISCUSSION 

We  performed  the  simulation  study  on  the  2NO^(NO)2 
system  at  low  P’s  and  T s  from  110  to  170  K,  and  for  nan¬ 
opore  widths  H=(0. 8, 1.0, 1.4, 1.7, 2.0)  nm.  The  experimen¬ 
tal  normal  boiling  point  of  the  bulk  system  is  121.4  K  and  at 
this  condition  the  system  is  almost  completely  monomeric  in 
the  vapor  phase  in  contrast  to  many  associating  fluids  such  as 
acetic  acid.  At  the  studied  P’s  and  T’s,  experiments  and 
RxMC  simulations  show  that  the  saturated-vapor  phase  con¬ 
tains  about  0.01-0.02  mole  fraction  of  dimers.  In  contrast  to 
the  saturated-vapor  phase,  the  dimerization  is  pronounced  in 
the  liquid  phase,  as  seen  from  Fig.  2.  Figure  2  shows  the 


FIG.  2.  Mole  fraction  of  (NO),  dimers  as  a  function  of  temperature  T  along 
the  saturation  line  for  the  bulk  liquid-phase  2NO=^ (NO),  system  from  the 
RxMC  simulations  and  from  the  experiments.  The  lines  are  drawn  as  a  guide 
to  the  eye. 


mole  fraction  of  liquid  (NO)2  along  the  experimental  satura¬ 
tion  line  (P  varies  from  0.2  to  43.3  bars),  as  obtained  from 
the  experiments  of  Smith  and  Johnston,  from  the  RxMC 
simulations  of  Turner  et  al.6  and  Johnson  et  al.,14  and  from 
our  RxMC  simulations.  Figure  2  shows  nearly  complete 
dimerization  at  low  T’s.  Further,  the  degree  of  dimerization 
decreases  with  increasing  T  and  at  P>170  K,  the  system 
becomes  monomeric.  We  also  see  from  Fig.  2  that  despite  the 
use  of  the  same  molecular  model,  our  RxMC  simulations 
predict  a  smaller  degree  of  dimerization  in  comparison  with 
the  RxMC  simulation  results  of  Turner  et  al.6  and  Johnson  et 
al . 14  This  is  a  consequence  of  the  differences  in  determining 
the  equilibrium  constant  K,  as  discussed  in  Sec.  Ill  C.  Since 
Turner  et  al.6  and  Johnson  et  al.14  adjusted  K  to  the  experi¬ 
mental  saturated-liquid  composition  at  110  and  116  K,  re¬ 
spectively,  their  RxMC  simulation  results  agree  quite  well 
with  the  experimental  data  measured  by  Smith  and 
Johnston  at  the  temperature  interval  from  1 10  to  120  K.  As 
previously  noted  by  Johnson  et  al.14  and  which  is  also  evi¬ 
dent  from  a  comparison  of  the  RxMC  simulation  results  of 
Turner  et  al.6  and  Johnson  et  al.14  in  Fig.  2,  the  simulation 
results  are  extremely  sensitive  to  changes  in  K.  We  are  not 
aware  of  any  experimental  data  on  the  saturated-liquid 
2NO^(NO)2  system  at  T>  120  K  to  shed  more  light  on 
differences  between  these  particular  RxMC  simulation  re¬ 
sults.  However,  we  argue  that  the  equilibrium  constant 
should  be  determined  unambiguously,  as  done  in  this  work, 
by  means  of  the  most  accurate  data  available  in  the  standard 
thermochemical  tables. 


A.  Effects  of  nanopore  size 

We  investigated  the  influence  of  the  nanopore  width  on 
reaction  conversion  by  performing  two-phase  reaction  simu¬ 
lations,  which  includes  a  bulk  phase  at  pressure  Pbuik 
=  0.16  bar  in  equilibrium  with  an  adsorbed  phase  in  a  model 
carbon  slit  nanopore.  Temperatures  were  varied  from 
120  to  160  K.  Due  to  the  low  value  of  Ebulk,  the  RxMC  simu- 
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(a) 


FIG.  3.  (a)  Mole  fraction  of  (NO)2  dimers  and  (b)  average  fluid  density  as 
a  function  of  temperature  T  in  the  model  carbon  slit  nanopore  of  various 
widths  H  at  a  bulk  pressure  of  0. 1 6  bar.  For  comparison,  the  simulation 
results  for  the  bulk  vapor  and  liquid  phases  are  also  included.  The  lines  are 
drawn  as  a  guide  to  the  eye. 

lations  showed  that  the  bulk  phase  behaves  as  an  ideal-gas 
system  and  is  predominantly  monomeric  with  about  a  0.01- 
0.02  mole  fraction  of  dimers. 

Figure  3  reports  the  yield  of  dimers  and  the  average  fluid 
density  in  the  pore  for  the  confined  NO  dimerization  reaction 
as  a  function  of  the  nanopore  width  and  temperature.  For 
comparison,  the  simulation  results  for  the  bulk  vapor  and 
liquid  phases  at  P= 0.1 6  bar  are  included.  Figure  3(a)  shows 
a  large  enhancement  of  the  dimer  formation  in  the  nanopores 
in  comparison  with  the  corresponding  bulk  vapor  phase.  The 
enhancement  becomes  greater  as  the  nanopore  width  is  re¬ 
duced.  For  the  nanopores  with  H 3=  1.7  nm,  the  influence  of 
further  increasing  H  becomes  small.  The  previous  studies  of 
Turner  et  al.6  and  Tripathi  and  Chapman18  resulted  in  similar 
findings. 

The  enhanced  conversion  is  a  result  of  two  phenomena 
occurring  simultaneously  in  the  nanopore  phase:  (i)  the  in¬ 
crease  in  the  average  fluid  density  in  the  pore  and  (ii)  the 
preferential  adsorption  of  the  (NO)2  dimer  on  the  nanopore 
surface.  The  enhanced  conversion  caused  by  the  higher  fluid 
density  in  the  pore  drives  the  reaction  equilibrium  towards 
the  formation  of  more  dimers,  as  seen  in  Fig.  3.  For  example. 


the  density  of  the  bulk  vapor  phase  at  7=  120  K  is  about 
0.16  X  10  4  mol/cm3  while  the  corresponding  average  fluid 
density  in  the  nanopores  with  H= 0.8  and  2.0  nm  is  approxi¬ 
mately  0. 14X1  0_1  and  0.87  X  10-2  mol/cm3,  respectively. 
On  the  other  hand,  the  average  fluid  density  in  the  nanopores 
is  substantially  lower  than  the  density  of  the  bulk  liquid 
phase  while  the  reaction  conversion  in  the  nanopores  with 
H= 0.8  and  1.0  nm  is  comparable  with  that  for  the  bulk  liq¬ 
uid  phase.  This  behavior  is  driven  by  the  preferential  adsorp¬ 
tion  of  the  dimer  in  the  pore  phase.  The  enhancement  is 
further  elucidated  in  plots  of  the  solid-fluid  interaction  po¬ 
tentials  Msf  in  Fig.  1,  in  profiles  of  the  NO  and  (NO)2  densi¬ 
ties  given  in  Fig.  4,  and  in  simulation  snapshots  shown  in 
Fig.  5.  As  seen  in  Fig.  1,  the  overlap  of  the  two  strongly 
adsorbing  wall  potentials  considerably  lowers  the  energy  in 
the  nanopore.  This  effect  is  more  pronounced  in  smaller  na¬ 
nopores,  where  values  of  usf  near  the  nanopore  center  drop 
well  below  those  for  the  larger  nanopores.  Hence  for  smaller 
nanopores,  the  fluid  molecules  near  the  nanopore  center  ex¬ 
perience  a  strong  interaction  with  both  walls.  Consequently, 
the  fluid  density  in  the  nanopore  center  is  high  in  these  small 
nanopores  [see,  e.g.,  the  fluid  density  profiles  and  simulation 
snapshot  for  H= 0.8  nm  in  Figs.  4(a)  and  5(a),  respectively], 
contributing  to  a  significantly  higher  dimer  yield.  The  pro¬ 
files  of  fluid  densities  in  the  smaller  nanopores  are  more 
dependent  on  the  variation  in  H  due  to  the  overlap  of  solid- 
fluid  interaction  potentials  near  the  nanopore  center,  as  seen 
for  H= 0.8,  1.0,  and  1.4  nm  in  Figs.  4(a)-4(c)  and  5(a)-5(c). 
Conversely  for  larger  nanopores,  the  fluid  molecules  near  the 
nanopore  center  experience  a  weak  interaction  with  the  walls 
and,  thus,  the  fluid  density  near  the  nanopore  center  is  small, 
as  seen  in  Figs.  4(d)  and  5(d)  for  the  nanopore  with  H 
=  2.0  nm.  It  should  be  also  mentioned  that  the  profiles  of 
fluid  densities  for  the  nanopore  with  H=  1 .7  nm  (not  shown 
here)  are  very  similar  to  those  for  the  nanopore  with  H 
=  2.0  nm.  Hence  for  these  larger  nanopores,  the  effect  of  the 
confinement  on  the  dimer  yield  is  less  strong  and  the  dimer 
yield  is  not  significantly  affected  by  the  variation  in  the  na¬ 
nopore  width. 

The  other  phenomena  occurring  in  the  nanopore  phase 
that  causes  an  enhanced  conversion  is  the  preferential  ad¬ 
sorption  of  the  (NO)2  dimer  on  the  nanopore  surface.  In 
other  words,  since  the  (NO)2  dimer  more  strongly  adsorbs 
inside  the  nanopore,  the  effect  seen  on  the  reaction  equilib¬ 
rium  appears  as  a  shift  to  a  higher  dimer  formation  inside  the 
nanopore.  Recall  that  the  interaction  energy  between  the 
(NO)2  dimer  and  the  nanopore  surface  is  approximately 
twice  as  large  as  the  NO  monomer.  However,  since  the  size 
of  the  (NO)2  dimer  is  approximately  twice  as  large  as  the 
NO  monomer,  the  relative  sizes  of  the  molecules  can  play  a 
countereffect.  The  extent  to  which  these  effects  balance  can 
be  seen  from  a  closer  examination  of  Figs.  4  and  5.  In  Figs. 
4(a)  and  5(a)  for  //=(). 8  nm,  it  is  clear  that  the  (NO)2  dimer 
strongly  adsorbs  in  the  center  of  the  nanopore  as  a  result  of 
the  influence  of  both  nanopore  surfaces.  At  this  nanopore 
size  there  is  little  remaining  available  space  in  the  nanopore 
for  the  NO  monomer  which  can  only  adsorb  on  the  nanopore 
surface.  As  the  nanopore  width  is  increased  to  H=  1.0  nm, 
the  profiles  of  fluid  densities  change  dramatically,  as  seen  in 
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FIG.  4.  Profiles  of  the  NO-  and  (NO)2-reduced  densities,  pa3,  in  the  model  carbon  slit  nanopore  with  widths  of  (a)  //-(). 8  nm,  (b)  //=].()  nm,  (c)  H 
=  1.4  nm,  and  (d)  H= 2.0  nm  at  a  temperature  of  125  K  and  at  a  bulk  pressure  of  0.16  bar;  z  is  the  distance  from  the  center  of  the  nanopore  which  is  located 
at  .-=().  Due  to  nanopore  symmetry,  only  half  of  the  fluid  density  profiles  are  shown. 


Figs.  4(b)  and  5(b).  The  larger  nanopore  width  allows  for  the 
formation  of  two  layers.  Now  both  species  form  layers  on  the 
nanopore  surface  with  the  more  strongly  attracted  (NO)2 
dimer  being  the  dominant  species.  From  Figs.  4(c)  and  5(c), 
we  can  see  that  as  the  nanopore  width  is  further  widened  to 
H=  1 .4  nm,  a  third  fluid  layer  emerges  at  the  center  of  the 
nanopore.  Flere  again  the  (NO)2  dimer  is  the  species  which 
more  strongly  adsorbs  shifting  the  reaction  equilibria  to  the 
right.  For  the  cases  where  H=  1.7  and  2.0  nm,  the  profiles  of 
fluid  densities  again  change  dramatically,  as  seen  from  Figs. 
4(d)  and  5(d)  for  H= 2.0  nm.  The  nanopore  is  now  so  wide 
that  fluid  particles  at  the  center  of  the  nanopore  do  not  expe¬ 
rience  the  attraction  of  both  surfaces,  thus  the  third  layer 
seen  in  Figs.  4(c)  and  5(c)  disappears.  The  profiles  of  fluid 


densities  in  Fig.  4(d)  have  well-defined  peaks  and  are  quite 
similar  for  both  species.  This  is  reflected  in  the  equilibrium 
concentrations  shown  in  Fig.  3(a). 

The  (NO)2  yield  in  the  nanopores  is  strongly  influenced 
by  the  temperature  (especially  for  //=1.4nm),  as  further 
seen  from  Fig.  3.  The  reaction  equilibrium  conversion  falls 
with  increasing  temperature,  as  observed  in  Fig.  3(a),  since 
the  NO  dimerization  reaction  is  exothermic  and  since  an  in¬ 
crease  in  T  is  accompanied  by  a  decrease  in  the  average  fluid 
density  in  the  pore,  as  shown  in  Fig.  3(b).  This  temperature 
effect  is  more  pronounced  than  that  seen  in  the  bulk  liquid 
phase. 

Finally,  Fig.  6  shows  the  one-dimensional  self-diffusion 
coefficients,  DN0  and  /J(NOh,  determined  from  NVT  MD 
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FIG.  5.  Snapshots  from  RxMC  simulations  of  the  NO/(NO)2  reactive  mixture  in  the  model  carbon  slit  nanopore  with  widths  of  (a)  7/=0.8  nm,  (b)  H 
=  1.0  nm,  (c)  H=\A  nm,  and  (d)  H= 2.0  nm  at  a  temperature  of  125  K  and  at  a  bulk  pressure  of  0.16  bar.  NO  and  (NO)2  molecules  are  shown  in  white  and 
red,  respectively.  The  snapshots  correspond  to  density  profiles  shown  in  Fig.  4. 
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FIG.  6.  One-dimensional  self-diffusion  coefficients  for  the  NO  monomers 
and  (NO)2  dimers,  Z>NO  and  Z)<NO*2,  in  the  model  carbon  slit  nanopore  as  a 
function  of  the  nanopore  width  H  at  a  temperature  of  125  K  and  at  a  bulk 
pressure  of  0.16  bar.  The  lines  are  drawn  as  a  guide  to  the  eye. 

simulations  of  nonreactive  mixtures  in  nanopores  of  various 
widths  at  7  =  1 25  K.  As  outlined  in  Sec.  II  B,  the  composi¬ 
tions  are  obtained  from  RxMC  simulations  and  held  constant 
during  the  NVT  MD  simulation.  We  see  from  Fig.  6  that  at 
77^  1.0  nm,  where  (NO)  2  is  the  abundant  component  of  the 
mixture  and  where  the  average  fluid  density  in  the  pore  in¬ 
creases  with  increasing  77,  values  of  Dsn  and  //N0)2  are 
approximately  the  same  and  change  only  slightly  with  in¬ 
creasing  77. 

As  the  nanopore  width  is  increased  to  1 .7  nm,  both  DN0 
and  Z)(n°)2  increase  where  the  effect  is  greater  for  NO  than 
for  (NO)2.  Here,  the  mole  fraction  of  (NO)2  as  well  as  the 
average  fluid  density  in  the  pore  decreases  with  increasing  77. 
Finally,  the  self-diffusion  coefficients  for  both  species  do  not 
change  significantly  between  77=1.7  and  2.0  nm  when  the 
system  becomes  mostly  monomeric.  The  overall  behavior 
exhibited  in  Fig.  6  is  indicative  of  the  shifting  reaction  equi¬ 
libria  with  respect  to  the  nanopore  width. 

B.  Impact  of  bulk  pressure  and  capillary  condensation 

We  studied  the  influence  of  the  bulk  pressure  and  capil¬ 
lary  condensation  on  the  NO  dimerization  reaction  equilib¬ 
rium  in  model  carbon  slit  nanopores  with  widths  of  1.4,  1.7, 
and  2.0  nm  at  T-  1 25  K.  Bulk  pressures,  Pbuik,  were  varied 
from  0.16  bar  up  to  the  experimental  vapor  pressure  of 
1.50  bar.  Capillary  condensation  is  a  first-order  transition 
that  leads  to  the  filling  of  nanopores  with  a  phase  character- 
ized  by  a  liquidlike  density."  We  found  that  the  nanopore 
with  77=  1.4  nm  is  too  small  to  exhibit  capillary  condensation 
while  the  nanopores  with  77=  1.7  and  2.0  nm  experience  cap¬ 
illary  condensation  around  Pbuik~0.4bar  and  TWk 
~  0.7  bar,  respectively.  These  results  are  discussed  in  more 
detail  below. 

1.  Model  carbon  slit  nanopore  with  a  width  of  1 .4  nm 

Figure  7  shows  the  mole  fraction  of  (NO)2  dimers,  the 
average  fluid  densities  in  the  pore,  and  the  one-dimensional 
self-diffusion  coefficients  as  functions  of  the  bulk  pressure 


FIG.  7.  (a)  Mole  fraction  of  (NO)2  dimers,  (b)  average  fluid  densities  in  the 
pore,  and  (c)  one-dimensional  self-diffusion  coefficients  for  the  NO  mono¬ 
mers  and  (NO)2  dimers,  Z)N0  and  /),  NO)2,  as  a  function  of  bulk  pressure  Pbulk 
in  the  model  carbon  slit  nanopore  with  a  width  of  1 .4  nm  at  a  temperature  of 
125  K.  The  symbols  are  simulation  results  of  this  work  and  the  lines  are 
drawn  as  a  guide  to  the  eye. 

for  the  confined  2NO^(NO)2  system.  We  see  from  Fig.  7(a) 
that  the  reaction  conversion  increases  with  increasing  Pbuik 
and  that  the  increase  is  rather  steep  between  0.16  and 
0.5  bar.  The  enhancement  of  the  reaction  conversion  is 
driven  by  an  increase  in  the  overall  average  fluid  density  in 
the  pore,  as  seen  from  Fig.  7(b).  The  enhancement  in  the 
(NO)2  yield  between  Pbulk=0.16  bar  and  Pbulk=1.50  bar  is 
about  20%,  while  the  corresponding  increase  in  the  overall 
average  fluid  density  in  the  pore  is  about  10%.  The  corre¬ 
sponding  7)n0  and  7)!XO)2  values  plotted  in  Fig.  7(c)  exhibit  a 
decrease  with  increasing  Pbulk,  i.e.,  with  an  increasing  overall 
average  fluid  density  in  the  pore.  The  dependence  of  Dh0  on 
P bulk  is  the  same  as  that  for  the  average  NO  density  in  the 
pore,  while  the  dependence  of  7) 1X0,2  on  Pbu]k  is  inverse  to 
that  of  the  average  (NO)2  density  in  the  pore  [cf.  Fig.  7(c) 
with  Fig.  7(b)].  The  values  of  l)K0  are  higher  than  those  of 
7><n°)2  despite  the  abundance  of  (NO)2  in  the  reactive  mix¬ 
ture  at  125  K  for  these  Pbuik’s.  Figure  7(a)  confirms  the  va- 
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FIG.  8.  Profiles  of  the  NO-  and  (NO)2-reduced  densities,  per3,  in  the  model 
carbon  slit  nanopore  with  a  width  of  1 .4  nm  at  a  temperature  of  125  K  and 
at  bulk  pressures  of  (a)  0.30,  (b)  0.75,  and  (c)  1.50  bar;  z  is  the  distance 
from  the  center  of  the  nanopore  which  is  located  at  c=0.  Due  to  nanopore 
symmetry,  only  half  of  the  fluid  density  profiles  are  shown. 


lidity  of  applying  Le  Chatelier’s  principle  to  explain  the  be¬ 
havior  of  the  NO  dimerization  reaction  equilibrium  in  slit 
nanopores:  the  equilibrium  of  a  reactive  system  in  which 
there  is  a  decrease  in  the  total  number  of  moles  shifts  toward 

products  when  such  a  system  experiences  an  increase  in 

23 

pressure. 

The  impact  of  Pbulk  on  the  structure  of  the  confined  re¬ 
active  mixture  of  NO/(NO)2  is  shown  in  Fig.  8  where  we 
plot  profiles  of  NO  and  (NO)2  densities  in  the  pore  for  dif¬ 
ferent  values  of  Fbulk.  We  see  from  Fig.  8  that  Pbulk  does  not 
affect  the  overall  character  of  both  profiles,  i.e.,  a  pro¬ 
nounced  peak  at  —0.37  nm  corresponding  to  the  location  of 
the  solid-fluid  potential  minimum  and  subsequent  peaks 
which  are  a  result  of  steric  and  packing  effects.  The  values  of 
P buik  only  influence  the  heights  of  these  peaks,  especially  for 


Pbuik  (bar) 


pbuik  (bar) 


FIG.  9.  (a)  Mole  fraction  of  (NO)2  dimers,  (b)  average  fluid  densities  in  the 
pore,  and  (c)  one-dimensional  self-diffusion  coefficients  for  the  NO  mono¬ 
mers  and  (NO)2  dimers,  DN0  and  D(N0^2,  as  a  function  of  bulk  pressure  Pbuik 
in  the  model  carbon  slit  nanopore  with  a  width  of  1 .7  nm  at  a  temperature  of 
125  K.  The  symbols  are  simulation  results  of  this  work  and  the  lines  are 
drawn  as  a  guide  to  the  eye.  The  adsorption  and  desorption  data  points  are 
shown  as  filled  and  open  symbols,  respectively. 


(NO)2,  as  seen,  e.g.,  from  a  comparison  of  Figs.  8(a)  and 
8(c).  It  is  evident  that  as  the  pressure  increases  the  adsorbed 
layers  become  more  defined,  but  only  slightly.  The  stronger 
adsorbing  (NO)2  dimer  layers  the  nanopore  surfaces  with  a 
less  defined  third  layer  residing  at  the  center  of  the  nanopore, 
while  the  NO  monomer  adsorbs  near  the  nanopore  surface. 


2.  Model  carbon  slit  nanopores  with  widths  of  1.7  and 
2.0  nm 

Figures  9  and  10  report  the  yield  of  dimers,  the  average 
fluid  densities,  and  the  diffusion  coefficients  Z5NO  and  Z)(N0* 2 
as  a  function  of  Fbulk  for  the  2NO=i  (NO)2  system  in  nan¬ 
opores  with  H-X.l  and  2.0  nm.  For  both  nanopore  widths, 
we  observe  a  steep  vertical  rise  in  the  adsorption/desorption 
isotherms  for  the  dimer  mole  fractions  [Figs.  9(a)  and  10(a)] 
and  for  the  average  fluid  densities  in  the  pore  [Figs.  9(b)  and 
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FIG.  10.  (a)  Mole  fraction  of  (NO)2  dimers,  (b)  average  fluid  densities  in 
the  pore,  and  (c)  one-dimensional  self-diffusion  coefficients  for  the  NO 
monomers  and  (NO)2  dimers,  Z>N0  and  Z)<NO*2,  as  a  function  of  bulk  pres¬ 
sure  Pjjuik  in  the  model  carbon  slit  nanopore  with  a  width  of  2.0  nm  at  a 
temperature  of  125  K.  The  symbols  are  simulation  results  of  this  work  and 
the  lines  are  drawn  as  a  guide  to  the  eye.  The  adsorption  and  desorption  data 
points  are  shown  as  filled  and  open  symbols,  respectively. 

10(b)]  as  a  result  of  capillary  condensation.  Furthermore, 
hysteresis,  widely  regarded  as  the  signature  of  capillary 
condensation,"  is  clearly  evident  in  these  isotherms.  The  bulk 
pressure  at  which  capillary  condensation  occurs  is  lower  in 
the  smaller  nanopore,  while  the  hysteresis  is  broader  in  the 
larger  nanopore.  A  dependence  of  DN0  and  Z)(NOb  on  Pbu|k 
[Figs.  9(c)  and  10(c)]  was  also  found.  Discontinuities  in 
these  isotherms  further  suggest  the  occurrence  of  capillary 
condensation. 

Closer  inspection  of  Figs.  9  and  10  reveals  some  very 
interesting  characteristics  of  the  reactive  system  caused  by 
capillary  condensation.  First,  the  reactants  and  products  are 
in  roughly  equal  proportion  prior  to  capillary  condensation 
for  both  nanopore  widths,  as  seen  from  Figs.  9(a)  and  10(a). 
However,  after  the  onset  of  capillary  condensation,  the  ratio 
of  dimers  to  monomers  changes  to  approximately  7:3  while 
reaching  a  ratio  of  8:2  at  the  highest  Pbulk  considered.  An¬ 
other  interesting  characteristic  is  that  the  increase  in  the 
overall  average  fluid  density  in  the  pore  as  a  result  of  capil¬ 


lary  condensation  is  very  nonuniformly  distributed  between 
the  monomers  and  the  dimers,  as  observed  in  Figs.  9(b)  and 
10(b).  Beyond  capillary  condensation,  the  dependence  of  the 
average  fluid  densities  in  the  pore  on  /Jbuik  is  very  similar  to 
that  found  in  the  H=1.4  nm  nanopore  [cf.  Figs.  9(b)  and 
10(b)  with  Fig.  7(b)].  Finally,  quite  unexpected  behavior  is 
the  dependence  of  DN0  on  Pbulk  prior  to  capillary  condensa¬ 
tion,  as  seen  in  Figs.  9(c)  and  10(c).  Here,  DN0  sharply  in¬ 
creases  with  increasing  Pbulk,  while  D(N0^2  decreases  with 
increasing  Pbulk.  Beyond  the  onset  of  capillary  condensation, 
the  dependence  of  Ds0  and  Z)(NOb  0n  Pbulk  resembles  the 
behavior  found  in  the  H=  1.4  nm  nanopore  [cf.  Figs.  9(c)  and 
10(c)  with  Fig.  7(c)], 

The  influence  of  capillary  condensation  on  the  structure 
of  the  confined  reactive  mixture  of  NO/(NO)2  is  elucidated 
in  Figs.  11  and  12.  Here,  we  display  the  profiles  of  the  NO 
and  (NO)2  densities  and  simulation  snapshots  in  the  nanop¬ 
ore  width  H=  2.0  nm  at  two  bulk  pressures  prior  to  capillary 
condensation,  T’buik=0-30  and  0.75  bar  [Figs.  11(a),  11(b), 
12(a),  and  12(b)],  and  at  two  bulk  pressures  after  the  onset  of 
capillary  condensation,  Fbuik=0.85  and  1.50  bar  [Figs.  11(c), 
11(d),  12(c),  and  12(d)].  It  should  be  noted  that  the  fluid 
density  profiles  in  the  nanopore  width  H=  1.7  nm  (not  shown 
here)  closely  resemble  fluid  density  profiles  in  the  H 
=  2.0  nm  nanopore.  Figures  11  and  12  clearly  show  the  ap¬ 
pearance  and  evolution  of  subsequent  dense  (condensation) 
layers  of  fluid  molecules  in  the  nanopore.  The  layer  posi¬ 
tioned  around  0.7  nm  from  the  nanopore  center  corresponds 
to  the  location  of  the  solid-fluid  potential  minimum.  The  ap¬ 
pearance  of  these  subsequent  dense  layers  is  manifested  by 
the  presence  of  additional  peaks  in  the  fluid  density  profiles. 
For  the  precondensation  case  Fbulk=0.30  bar  [Figs.  11(a)  and 
12(a)],  the  onset  of  another  layer  of  NO  molecules  between 
0.3  and  0.4  nm  from  the  nanopore  center  is  evident.  Then 
just  prior  to  capillary  condensation,  i.e.,  at  Fbulk=0.75  bar, 
this  additional  NO  layer  becomes  further  pronounced  while 
the  onset  of  an  additional  layer  of  (NO)2  molecules  posi¬ 
tioned  also  between  0.3  and  0.5  nm  from  the  nanopore  center 
can  be  seen  in  Figs.  11(b)  and  12(b).  Finally,  after  capillary 
condensation,  these  additional  NO  and  (NO)2  layers  further 
evolve  along  with  the  appearance  of  layers  of  NO  and  (NO)2 
molecules  in  the  nanopore  center,  as  seen  from  Figs.  11(c), 
11(d),  12(c),  and  12(d). 

V.  CONCLUSIONS 

The  influence  of  confinement  on  chemical  reaction  equi¬ 
librium  in  nanoporous  materials  was  studied  in  detail  by 
means  of  reaction  ensemble  Monte  Carlo  and  molecular- 
dynamics  simulations  for  the  nitric  oxide  dimerization  reac¬ 
tion  in  model  carbon  slit  nanopores.  The  effects  of  tempera¬ 
ture,  slit  width,  bulk  pressure,  and  capillary  condensation  on 
the  reaction  conversion,  fluid  structure,  and  one-dimensional 
self-diffusion  coefficients  have  been  reported  and  discussed 
thoroughly. 

A  large  increase  in  the  dimerization  in  the  nanopore 
phase  was  found  with  respect  to  the  dimerization  in  the  cor¬ 
responding  bulk  phase.  The  dimerization  increased  by  a  fac¬ 
tor  of  more  than  80  for  the  smaller  nanopores  at  the  lower 
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FIG.  11.  Profiles  of  the  NO-  and  (NO)2-reduced  densities,  po 3,  in  the  model  carbon  slit  nanopore  with  a  width  of  2.0  nm  at  a  temperature  of  125  K  and  at 
bulk  pressures  of  (a)  0.30  bar,  (b)  0.75  bar  (prior  to  capillary  condensation),  (c)  0.85  bar  (after  capillary  condensation),  and  (d)  1.50  bars;  z  is  the  distance 
from  the  center  of  the  nanopore  which  is  located  at  z=0.  Due  to  nanopore  symmetry,  only  half  of  the  fluid  density  profiles  are  shown. 


temperatures  and  low  pressure.  The  enhanced  dimerization  is 
due  to  the  combined  effects  of  the  increased  fluid  density  in 
the  nanopore  phase  and  the  preferential  adsorption  of  the 
(NO)2  dimer  in  the  nanopore.  Analogous  to  the  bulk  phase 
behavior,  the  dimerization  decreases  with  increasing  tem¬ 
perature.  The  impact  of  the  bulk  phase  pressure  on  the  reac¬ 
tion  conversion  is  moderate,  unless  capillary  condensation 


occurs  in  the  nanopores,  in  which  case  then  is  quite  dramatic. 
The  capillary  condensation  changes  the  ratio  of  dimers  to 
monomers  from  roughly  1:1  prior  to  capillary  condensation 
to  a  ratio  of  7:3  after  capillary  condensation.  This  enhanced 
dimerization  is  caused  by  an  increase  in  the  overall  fluid 
density  in  the  nanopores.  However,  this  increase  is  not  dis¬ 
tributed  uniformly  between  monomers  and  dimers.  There  is  a 
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FIG.  12.  Snapshots  from  RxMC  simulations  of  the  NO/(NO)2  reactive 
mixture  in  the  model  carbon  slit  nanopore  with  a  width  of  2.0  nm  at  a 
temperature  of  125  K  and  at  bulk  pressures  of  (a)  0.30  bar,  (b)  0.75  bar 
(prior  to  capillary  condensation),  (c)  0.85  bar  (after  capillary  condensation), 
and  (d)  1.50  bar.  NO  and  (NO)2  molecules  are  shown  in  white  and  red, 
respectively.  The  snapshots  correspond  to  density  profiles  shown  in  Fig.  1 1 . 


significantly  larger  increase  in  the  average  dimer  density  in 
the  pore  than  in  the  average  monomer  density  in  the  pore 
upon  capillary  condensation. 

The  magnetic  susceptibility  measurements  of  Kaneko  et 
al.9  and  Nishi  et  al.10  have  found  a  much  larger  increase  in 
dimerization  due  to  confinement  than  the  simulation  studies.6 
Although  the  experiments  and  simulations  in  this  work  are 
performed  at  different  thermodynamic  conditions,  some 
qualitative  comparisons  can  be  made.  First,  discrepancies 
may  be  due  to  experimental  challenges  in  distinguishing  be¬ 
tween  equilibrium  states  and  long-lived  metastable  states  in 
the  nanopore.  Second,  other  molecular  species  formed  within 
the  nanopore  may  be  incorrectly  interpreted  as  (NO)2  dimers 
including  complexes  formed  between  the  NO  monomer  and 
defects  on  the  nanopore  surface.  While  the  intermolecular 
potential  models  used  in  the  simulations  are  idealized, 
changes  in  the  carbon-adsorbate  interactions  would  need  to 
be  unrealistically  large  to  match  the  experimental  findings. 
For  example,  Turner  estimated  that  the  binding  energy  of  the 
dimer  would  have  to  increase  by  a  factor  of  6  in  order  to 
approach  the  elevated  conversions  measured  in  the 
experiments.43  An  increase  in  the  binding  energy  of  this 
magnitude  is  unlikely.  Moreover,  the  excellent  quantitative 
agreement  between  simulation  and  DFT  calculations13  pro¬ 
vides  further  justification  for  reconsidering  the  experimental 
results. 
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